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ABSTRACT 

High redshift quasars can be used to trace the early growth of massive galaxies and may 
be triggered by galaxy-galaxy interactions. We present MUSE science verification data 
on one such interacting system consisting of the well-studied z = 3.2 PKS1614+051 
quasar, its AGN companion galaxy and bridge of material radiating in Lya between 
the quasar and its companion. We find a total of four companion galaxies (at least two 
galaxies are new discoveries), three of which reside within the likely virial radius of 
the quasar host, suggesting that the system will evolve into a massive elliptical galaxy 
by the present day. The MUSE data are of sufficient quality to split the extended Lya 
emission line into narrow velocity channels. In these the gas can be seen extending 
towards each of the three neighbouring galaxies suggesting that the emission-line gas 
originates in a gravitational interaction between the galaxies and the quasar host. 
The photoionization source of this gas is less clear but is probably dominated by 
the two AGN. The quasar’s Lya emission spectrum is double-peaked, likely due to 
absorbing neutral material at the quasar’s systemic redshift with a low column density 
as no damping wings are present. The spectral profiles of the AGN and bridge’s Lya 
emission are also consistent with absorption at the same redshift indicating this neutral 
material may extend over > 50 kpc. The fact that the neutral material is seen in the 
line of sight to the quasar and transverse to it, and the fact that we see the quasar and 
it also illuminates the emission-line bridge, suggests the quasar radiates isotropically 
and any obscuring torus is small. These results demonstrate the power of MUSE for 
investigating the dynamics of interacting systems at high redshift. 

Key words: galaxies: interactions - quasars: individual: PKS1614+051 - galaxies: 
high-redshift - techniques: imaging spectroscopy 


1 INTRODUCTION 

The growth of galaxies and galaxy clusters can be probed 
using high redshift qu asars to pinpoint regions of early 
clu ster f ormat ion (e.g . iTurnerl Il99ll : ISpringel et al.l 120051 ; 
IDi Matteo et alj 120081 1. The rapid buildup in mass of the 
central supermassive black holes and host galaxies of quasars 
implies that they are hosted by massive dark matter haloes 
- themselves the likely sites of early galaxy clustering. Thus 
distant quasars can be used to explore the early inter¬ 
actions that may trigger quasar activity and lead to the 
growth of massive galaxies that eventually exist at the cen¬ 
tres of present-day clustering. Observational evidence of 
galaxy clustering around at least a subset of high redshift 
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quasars is now becoming apparent (e.g. Toshikaw a et al.l 
120121 ; iHusband et al.ll2013 : iMcGreer et al.ll2()14fh 


One of the earliest identified systems that showed evi¬ 
dence of interaction and clustering at high redshift is that 
of PKS1614+051, consisting of a z = 3.2 radio-loud quasar 
and companion galaxy connected by an extended bridge 
of ionized gas . The emission-line companion galaxy is it¬ 
self a n AGN llDiorgovski et all 1 19871 : iBremer fe Johnstonel 
Il995lf an d the bri dge of gas is radiating strongly in Lya 
( Hu fe Cowielll987l l. While not unique in having clearly ex¬ 
tended Lya, the combination of quasar, AGN and extended 
Lya emission is a rare opportunity to explore the connec¬ 
tion between a hig h redshift quasar and i ts environment. 
First discovered by iDiorgovski et alj (1 1985h . the brightness 
of the companion galaxy (significantly brighter than typical 
Lyman break galaxies and Lya emitters at 2 ~ 3) and the 
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Lya bridge allowed the system to be studied with 4-m tele¬ 
scopes in the 1980s and 19 90s. Given the companion galaxy, 
some authors fe.g. lHu fc Cowidll987 l suggest it should lie in 
an overdense (or protocluster) environment that will evolve 
into a galaxy cluster by the present day. 

With the advent of MUSE (Multi-Object Spectroscopic 
Explorer; iBacon et al . m3) - a 1 x 1 arcmin 2 field inte¬ 
gral field unit (IFU) on the VLT (Very Large Telescope) - a 
more detailed study of this region is now possible. The large 
IFU enables continuum and emission line morphologies to 
be acquired as well as information about the dynamics and 
ionization state of the gas and galaxies, without which in¬ 
terpretation of the complex interaction of quasars with their 
environments is highly challenging. Often the close environ¬ 
ment of a quasar is studied using rare, spatially coincident 
background quasars to probe individual sight lines through 
the halo containing the quasar. Here we show that we can 
probe numerous lines of sight around a quasar using the IFU 
capabilities of MUSE and the Lya halo of the quasar. Al¬ 
though already well-studied, the MUSE science verification 
data presented in this paper provides new information on 
the nature of this system. 

A ACDM cosmology with Hp = 70 km s -1 Mpc -1 , 

= 0.29 and Ha = 0.71 dBennett et al.l 120141 ') is used 
throughout. 


2 OBSERVATIONS 

The data were taken as part of Science Verification with 
VLT/MUSE between 18 and 20th August 2014. The lxl 
arcmin 2 field of view was first centered 15 arcsec west of 
PKS1614+051 and then 15 arcsec to the east. Four 880s ex¬ 
posures, rotated 90 deg between each, were taken at each 
position, giving a total exposure time of 1.96 h and spectro¬ 
scopic coverage of all galaxies within 1.28 x 1 arcmin 2 over a 
wavelength range of 4750-9350 A . This arrangement meant 
that the central ~30 arcsec wide North-South strip of the 
eventual image accrued twice the exposure time of the ~15 
arcsec wide eastern and western strips. 

The MUSE Data Reduction Software (version 0.18.5) 
was used for the first reduction of the data - bias subtract¬ 
ing, flat fielding and calibrating the data, and aligning and 
constructing three dimensional data cubes. An IDL script 
was written to achieve a better sky subtraction than the 
early MUSE pipeline by fitting and subtracting a 2nd order 
polynomial to each spatial row (along the long axis of the 
IFUs) of each wavelength slice of the resulting data cubes 
for each of the eight exposure in order to remove the sky 
contribution. The resulting eight data cubes (one for each 
of the four rotations at two separate pointings) were mean- 
combined with cosmic ray rejection (sigma clipping) and 
the median of the entire frame at each wavelength bin sub¬ 
tracted to remove any remaining residual from the sky sub¬ 
traction. Images were constructed from the data cubes by 
collapsing the spectral axis over defined wavelength ranges. 
A further median subtraction was done to each row of each 
individual IFU for every exposure and then these were com¬ 
bined with a bad pixel mask and sigma clipping to produce 
the final white-light image (see Fig. 0. 

The RMS noise in the final data cube is 1 x 10~ 19 erg 
s' 1 cm 2 per 1.25 A wavelength bin, at 5150 A and in the 


center of the data cube. The RMS noise is a factor of %/2 
worse in the two 15 arcsec wide strips on the east and west 
of the image where the total exposure time is half that of the 
middle and the noise is up to 50 per cent lower in the mid¬ 
dle of the wavelength range around 7000 A. The 2<r depth, 
within 2 arcsec diameter apertures, of the white light, full 
range (4750-9350 A) stacked image is 26.4 AB. The FWHM 
velocity resolution is measured from the sky lines to be 2.7 
A or equivalently 150 km s -1 at 2 = 3.215 and the spatial 
FWHM of the final white light image is 0.9 arcsec. The pixel 
scale of MUSE is 0.2 arcsec (pixel) -1 . 

The final white light image is sho wn in Fig. [T] We 
use SExtractor dBertin fc ArnoutsHl996l) followed by a vi¬ 
sual check to identify 119 objects in our white-light image, 
of which 60 are brighter than R = 25 AB (see Table |T|). 
We also checked the entire data cube for emission lines but 
only one additional object was found. This is unsurprising as 
any emission line only object would have an extremely high 
equivalent width and likely a very red or very blue spectrum 
in order for the continuum to remain undetected. The only 
other object with undetected continuum is the Lya bridge 
indicating it has a very large rest-frame equivalent width 
(>60 A using the 2<r limit on the spectrum in the wavelength 
range 5090-5190 A near the Lya emission or >180 A using 
the 2a limit on the full range white light image shown in 
Fig- 0 that is unlikely to be powered by embedded star for¬ 
mation in the bridge (the max imum equivalent width from 
star forming regions is 240 A: iHennawi fc Prochaskall2013l : 
ICharlot fc Falll 119931 '). Analysis of both emission lines, ab¬ 
sorption lines and comparison with Sloan Digital Sky Sur¬ 
vey (SDSS) galaxy and stellar templates results in redshifts 
for 34 galaxies (Table QJ, of which four are at the redshift 
of the quasar (within Az = 0.003). We are sensitive to Lya 
flux densities of greater than 2 x 10 -19 erg s _1 cm -2 A -1 
at the redshift of the quasar, determined by collapsing a 
100 A wide section of the data cube around 5120 A over 
the wavelength axis. In addition, we found no convincing 
Lyman break galaxies (LBGs) without Lya emission when 
we compared stacked images above and below potential Lya 
breaks at 20 A intervals over the full wavelength range. 


3 RESULTS 

The capabilities of MUSE are clearly demonstrated in this 
data - with just two hours of observations we can identify 
galaxies at the same redshift as the quasar across the 1 
arcmin 2 field of view and probe in detail the region close 
to the quasar containing PKS1614, its AGN companion, 
the emission line bridge and any other potential compan¬ 
ion galaxies. 


3.1 Galaxy Environment of PKS1614 

The PKS1614+051 quasar has at least four companion 
galaxies: the previously known narrow line AGN A, galaxy 
E - hitherto identified as a candidate z = 3.2 L ya emitter 
(LAE) based o n narrow-band imaging dHu fc Cowiel 1 19871 : 
iHu et al.| fl99l]) - and two new LAEs/LBGs: G and BH. For 
these new ly identified sources , we extend the galaxy labeling 
scheme of lHu fc Cowiel dl987f ) but we do not find galaxies at 
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Figure 1. A white light image (4750-9350 A) of the MUSE field of view highlighting the galaxies whose redshift could be determined as 
well as the stars that could be identified. Galaxies that are at the same redshift as the central quasar (Q) are highli ghted in green (A, 
E, G and BH), whilst those at lower redshift are in red. The labeling for the central objects is the same as in Hu Sz Cowi c lll987T ) with 
Q marking the quasar. All the details for these objects are in Table [T] 


the positions of their C and D objects. The spectra of these 
companion galaxies are shown in more detail in Fig. [2]&[3] 

Three of the companion galaxies (A, E and G) are 
clo se enoug h to be int eracting with the quasar as defined 
by lEllison et al.l (l2008lf and are within the expected virial 
radius of the quasar’s host dark matter halo suggesting they 
will eventually merge with the quasar host to form a massive 
elliptical galaxy (perhaps at the centre of a group or clus¬ 
ter) by the present day. Although the spectrum of E may be 
affected by scattered light from the quasar, its broad-band 
photometry is consistent with a Lyman break and its Lya 
line is significantly narrower than the quasar’s, suggesting 
it is a separate a = 3.2 galaxy. E is unlikel y to be driven 
solely by fluorescence from t he quasar ('e.g. lCantalupo et al.l 
2005:; iKollmeier et al.l 120100 as the presence of continuum 
emission suggests star formation can create the Lya emis¬ 
sion (the rest-frame equivalent width is only ~15 A). The 
companion galaxies are all cleanly separated from the quasar 
demonstrating the excellent spatial sampling of the MUSE 
data. 


In total, we confirm five galaxies including the quasar 
host within an area of 20 by 25 arcsec by Az = 0.003 or a 
maximum co-moving volume of 1.4 co-moving Mpc 3 , ignor¬ 
ing peculiar velocities i.e. assuming the difference in redshift 
represents the Hubble expansion. 


3.2 Gaseous Environment 

We have obtained a detailed, high signal-to-noise spectrum 
of the quasar, companion A and, for the first time, the emis¬ 
sion line bridge. Fig. [4] shows the distribution of gas in the 
central region around the quasar. The emission line bridge 
extends over 50 physical kpc from A to the quasar and 
slightly beyond the quasar in the other direction (Fig. (4jr), 
with a width of around 10 kpc. Fig. [5}; shows the red-most 
and blue-most (5 A wide or 300 km s -1 ) components of the 
Lya line, showing how the high surface brightness emission 
connects the quasar to companion A and also extends to¬ 
wards the two other nearby galaxies (E and G), denoted by 
crosses. Fig. |4jl shows the dispersion of the Lya line as a 
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ID 

2 peak 

Zsys 

z fit 

Lya flux 

/10~20 er g s -l cm -2 

Lya Luminosity 
/10 43 erg s -1 

Lya Vel. Dispersion 
/km s —1 

R Abs. Mag. 
/mag 

Q 

3.215 

3.212 

3.215 

616720 

58.8 

920 

-28.1 

A 

3.215 

3.214 

3.215 

33260 

3.18 

400 

-23.0 

E 

3.215 

na 

3.215 

3960 

0.035 

260 

-23.2 

G 

3.215 

na 

3.211 

1080 

0.103 

190 

-22.7 

BH 

3.217 

na 

3.217 

773 

0.074 

130 

-22.9 

Bridge 

3.215 

na 

3.214 

10940 

1.989 

250 

-21.0 


Table 2. The properties of the companion galaxies relative to the quasar. The Lya flux is determined from the area under a Gaussian 
fit to the line with the continuum removed in a 2 arcsec aperture. The properties for the bridge are for the entire bridge structure once 
the contribution from the quasar and companion A are removed, except for the velocity dispersion which is the value for the centre of 
the bridge. 


ID 

RA 

Dec 

Redshift 

Classification 

Q 

16:16:37.56 

04:59:32.1 

3.212 

Quasar 

A 

16:16:37.83 

04:59:37.1 

3.214 

LAE/NLAGN 1 

B 

16:16:37.52 

04:59:22.9 

0.507 

Spiral galaxy 

E 

16:16:37.72 

04:59:33.8 

3.215 

LAE 

F 

16:16:37.97 

04:59:24.0 

0.67 

Galaxy 

G 

16:16:37.59 

04:59:36.2 

3.215 

LAE 

H 

16:16:37.34 

04:59:25.4 

0.000 

M star 

I 

16:16:38.49 

04:59:20.5 

0.143 

Galaxy 

J1 

16:16:36.69 

04:59:23.3 

0.53 

Early-type galaxy 

J2 

16:16:36.77 

04:59:21.3 

0.615 

Luminous red galaxy 

L 

16:16:38.11 

04:59:12.8 

3.098 

LBG/LAE 

N 

16:16:36.77 

04:59:12.3 

1.089 

Galaxy or 
Interacting galaxies 

O 

16:16:36.58 

04:59:13.1 

0.000 

M Star 

PI 

16:16:36.66 

04:59:09.5 

0.608 

Galaxy 

P2 

16:16:36.65 

04:59:08.1 

0.608 

Galaxy 

R 

16:16:39.20 

04:59:08.7 

0.000 

Star 

S 

16:16:39.52 

04:59:06.7 

0.402 

Galaxy 

T 

16:16:39.87 

04:59:23.3 

0.905 

Galaxy 

V 

16:16:40.11 

04:59:38.3 

0.71 

Galaxy 

W 

16:16:40.08 

04:59:35.9 

0.16 

Galaxy 

Y 

16:16:39.17 

04:59:53.3 

1.470 

Galaxy 

Z 

16:16:39.12 

04:59:55.9 

0.000 

M star 

AF 

16:16:36.54 

04:59:58.1 

0.000 

M star 

AG 

16:16:35.94 

05:00:01.3 

1.213 

Galaxy 

AH 

16:16:36.03 

04:59:55.7 

0.833 

Galaxy 

AI 

16:16:36.40 

04:59:45.1 

0.000 

Star 

AR 

16:16:35.80 

04:59:20.3 

0.629 

Galaxy 

AT 

16:16:35.11 

04:59:12.1 

0.549 

Spiral galaxy 

AV 

16:16:38.34 

04:59:04.5 

0.000 

Star 

AW 

16:16:37.35 

05:00:03.3 

0.000 

M star 

BD 

16:16:39.44 

04:59:16.9 

0.143 

Galaxy 

BH 

16:16:37.13 

04:59:54.5 

3.217 

LAE 

BM 

16:16:39.31 

04:59:20.7 

0.853 

Galaxy 

BN 

16:16:37.05 

04:59:35.3 

1.050 

Galaxy 


Table 1. A summary of the redshifts obtained. The ID is the 
same as in Fig. IT1 and is an extension of the labeling scheme of 
iHu &: Cowid d 1987l'> ■ W e do not find galaxies at the positions of 
their C and D objects, feremer fe Johnstonel lll995l i 

function of position, found from fitting a Gaussian to the 
Lya line in each spatial pixel. The Lya emission from the 
bridge has a relatively low dispersion, as does the Lya emis¬ 
sion from companions E and G, whilst the emission from the 
quasar and companion A have larger velocity dispersions as 
expected for active galaxies (see Tabic [2]). 

The spectrum of the quasar’s Lya emission is strongly 
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Figure 2. The extracted ID spectrum of companion A with the 
sky lines plotted underneath. 


double-peaked with the central gap having a rest-frame 
equivalent width of 1.1 A and the flux declining to almost 
zero at the central wavelength of the line. If this is due to 
an absorption system the velocity width, b, is 180±15 km 
s' 1 and this implies a low column density of (1.9 ± 0.4) x 
10 13 cm -2 using a curve of growth analysis and assuming 
a covering fraction of 100 per cent, as expected for an un¬ 
saturated absorption profile without damping wings. This 
potential absorbing material is at the systemic redshift of 
the quasar (as derived from other lines e.g. Hell, CIV) sug¬ 
gesting that the emission arises from neutral material in or 
around the quasar. The spectra of A and the bridge are also 
consistent with being affected by the same absorbing ma¬ 
terial. Their spectral profiles display a similar sharp cutoff 
on the blue-ward edge as the quasar’s Lya profile (Fig. 0, 
although the underlying Lya emission lines in A and the 
bridge are narrower and appear to peak at slightly longer 
wavelengths than the quasar Lya emission line. 

We see no other emission lines arising from the bridge, 
placing a constraint on the CIV/Lya and Hell/Lya flux ra¬ 
tios of < 0.03 and < 0.06 (3cr), significantly below the ratios 
for the quasar and companion A where CIV and Hell was 
detected. This implies a lower ionization parameter for the 
bridge than that typically seen from an AGN (see discussion 
section below). 
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Figure 4. a. White light image of the centre of the MUSE field of view (20 by 18 arcsec) clearly showing galaxies E and G. b. The same 
area in the narrow band (5120-5135 A) encompassing Lya. Some of the Lya emission from the quasar and A has also been subtracted 
using a PSF model in order to show the extended emission more clearly. The contours trace the Lya showing that both the quasar and 
A are clearly extended in Lya and the emission extends beyond A to the North. The crosses mark the positions of galaxies E and G. The 
surface brightness limit is ~ 4 X Hr 17 erg s 1 cm -2 arcsec -2 and the contours indicate surface brightness of (4, 5, 9, 11 and 20) X10" 1 ' 
erg s -1 cm -2 arcsec -2 . c. The Lya line is split into the blue-most (5120-5125 A) and red-most (5130-5135 A) velocity components. 
In these narrower components the gas is seen to extended between the quasar and A and also towards galaxies E and G (marked by 
crosses), d. The velocity dispersion of the Lya line from the same Gaussian fits to each spatial pixel. This shows that the quasar has a 
velocity dispersion of 920 km s —1 whilst companion A has a velocity dispersion of 400 km s —1 . The high values are what is expected 
for active galaxies. Companion E and the material in the bridge have much lower velocity dispersions, around the velocity resolution of 
MUSE. 


4 DISCUSSION 

4.1 Protocluster or Protogalaxy Environment 

PKS1614+051 lies in a galaxy overdensity of at least five 
galaxies (including itself) within a volume of no more than 
1.4 co-moving Mpc 3 , demonstrating that MUSE is clearly 
capable of identifying clustering and potential protoclusters 
at z ~ 3. At this redshift, MUSE allows for little spectral 
coverage below the 1216 A Lya break. Given the relative 


weakness of this break at 2 = 3 when compared to higher 
redshifts, and the low signal to noise obtained on the contin¬ 
uum of faint sources at short wavelengths, it will be easier 
to use MUSE to search for Lyman break galaxies at higher 
redshifts than that probed here. 

The majority (3/4) of the companion galaxies are within 
a 50 kpc of the quasar, a nd hence well within it s expected 
virial radius (~ 280 koc: ICantahmo et al.l f2014). suggest¬ 
ing that they will merge with the quasar to form a massive 
elliptical galaxy by 2 = 0. If this is the progenitor of the 
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Figure 3. The extracted ID spectrum of companion E, G and 
BH with the sky lines plotted underneath. The dashed line shows 
the median of the spectrum either side of the expected break. 
Areas of the spectrum affected by skylines have been set to zero. 



Wavelength /Ang 

Figure 5. A zoom in on the extracted ID spectrum of the quasar, 
companion A and the bridge showing the Lya emission profile. 
The flux of the bridge and quasar have been scaled to that of 
companion A. There is strong absorption at 5120 A that is Lya 
at the quasar’s systemic redshift. All three objects show a sharp 
drop at this wavelength. The bridge appears to recover flux short- 
ward of the drop, as may companion A but note that the signal 
to noise on the spectrum of A is much higher than the signal to 
noise of the bridge. 


massive galaxies seen at z ~ 1.5 then the merger needs to 
be completed within ~ 2.5 Gyr. 

At this redshift we expect protoclust ers to extend over 
radii of 2-5 arcmin dChiang et al.l l2013t) . As MUSE only 
probes 1 arcmin 2 we can draw limited conclusions about 
the wider protocluster environment of PKS1614. Neverthe¬ 
less, with one companion galaxy outside its virial radius and 
three within, the quasar is highly likely to reside in an over- 
dense environment likely to evolve into a group or cluster. 
We note that with another eight similar pointings - in total 
18 hours of MUSE time - the environment of the system 
could be probed in detail on the scale of an ACS or WFC3 
HST image, better matched to the expected scale of over- 
dense regions. 


4.2 Gaseous Interactions 


High redshift radio loud quasars and galaxies are of¬ 
ten observed to have extended L ya haloes su rrou nding 
them, spanning lOs-lOOs of kpc (e.g. Heckman et al.lll991bl : 
IChristensen et ahll2006l : IVillar-Martm et al. 2007h with line 


widths of 100-1500 km s -1 I Heckman et al.l Il991al h The 
photoionization of the gas is often attributed to radi¬ 
ation escaping along the radio axis of the galaxy, al¬ 
though fainter extended Lya haloes are seen around com¬ 
pact radio galaxies (e.g. Ivan Oiik et al.| 1997 ) and a round 
some radio-quiet quasars (e.g. Christensen et all 120061 : 
iHerenz et al.ll2015n . Many of the extended emission line re¬ 
gions around radio galaxies (~ 60 per cent) and in par¬ 
ticular those around the smallest (< 50 kpc sized) radio 
sourc es such a s PKS1614+051, have associated HI absorp¬ 
tion (Ivan Oiik et al. 1997 )- Indeed ~ 60% of the 74 z ~ 2 
quasars studied bv~ IProchaska et al.l J2013bfo were found to 
have nearby (< 200 kpc) neutral material absorbing pho¬ 
tons from a spatially coincident background quasar. Assum¬ 
ing the quasar is at the centre of a similar large-scale diffuse, 
high filling factor neutral halo, the picture we favour to ex¬ 
plain the observed distribution of gas around PKS1614+051 
is as follows: Embedded within this diffuse halo are also the 
companion galaxies A, E and G and the bridge of denser, 
lower filling factor emission-line gas. This gas phase orig¬ 
inates from material extracted from the companions and 
quasar host through mutual gravitational interactions and 
is ionized by AGN and/or stellar emission. Below we discuss 
the evidence in favour of this picture. 

The gas in the bridge is aligned between the quasar and 
companion A, and may extend slightly beyond A and the 
quasar along the same axis. Fig. St shows the distribution 
of the red-most and blue-most component of the Lya bridge 
gas. There is a comparatively narrow stream of gas connect¬ 
ing the quasar to companion A, with a wavelength similar 
to the Lya peak wavelength of A, suggesting the two galax¬ 
ies are gravitationally interacting. The bridge also seems to 
extend towards galaxies E and G (marked by crosses), giv¬ 
ing it its broad appearance in the full-width Lya image and 
indicating that these galaxies may also be interacting with 
the bridge and/or quasar. 

Unlike many powerful high redshift radio sources, in 
this case there is no obvious spatial relationship between 
the radio emission and the extended emission-line region. 
Radio observations of the quasar indicate that it is a gi¬ 
gahertz peaked source with the bulk of the radio emission 
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arising from a milli-arcsec scale region. The source’s radio 
jet originates from the quasar in a south-westerly direction 
before turning throug h 90 degrees to the no r th-west again 
on su b arcsec scales dDiorgovski et all Il987l : lOrienti et al.l 
2006) and hence the extended emission line region is prob¬ 
ably aligned with the counter-jet axis, albeit on scales two 
or more orders of magnitude larger than that of the radio 
emission. However, the lack of similarly extended emission 
in the south-westerly direction suggest the radio-jet is not 
the primary cause of the emission-line bridge. 

The only emission line detected from the gas in the 
bridge is Lya, placing strong constraints on the ratio of 
CIV/Lya and Hell/Lya. The lack of Hell suggests a rel¬ 
atively low ionization parameter, unlike that seen in gas 
coincident with jets in radio galaxies, whereas the lack of 
CIV emission suggests the gas is not hot enough to ion¬ 
ize CIV. Although gas near AGN typically has a high 
ionization parameter, low ionization parameter gas may 
arise at greater distances from the AGN where the pho¬ 
ton field is more diffuse, particularly if the ionizing pho- 
tons have already passed t hrough a cloudy medium (e.g. 
lArrigoni Battaia et al.ll2015f ). The majority of ionizing pho¬ 
tons exciting the bridge emission probably emanate from 
the quasar and the AGN within companion A as the Lya 
surface brightness rises towards both (Fig. [6]), although the 
exact form of the surface brightness profile will depend on 
the optical thickness and density distribution of the gas in 
the bridge. The large line width of Lya suggests the bridge 
emission is not due to gravitational cooling, although grav- 
itational coo ling would produce little or no Hell and CIV 
(lYang et al.ll200a ) , and the Lya emission is unlikely to arise 
from shocks as this would result in a hi gher fraction of 
CIV c ompared to Lya than that seen here ( Heckman et al.l 
Il991al h Th e line ratios are cons istent with photoionization 
from stars dPrescott et al.ll2009l !. although the Lya equiva¬ 
lent width limit of the bridge is too high for the ionization 
source to be solely stellar (the equivalent width limit only 
allows for ionization from massive, young O stars and such 
a strong star burst would be a very short lived phenomenon 
making this an unlikely scenario). However, a more typical 
stellar population could provide some fraction of the ionizing 
photons along with ionizing photons from the two AGNs. Al¬ 
ternatively, the observed line ratios could just be an effect of 
the resonant properties of Lya whereby Lya will resonantly 
scatter out to large r distances from the quasar than non¬ 
resonant Hell (lArrigoni Battaia et al.l 120141 1 . However, this 
may not be a dominant effect as the photons that resonantly 


scatter are likely to escape the system on much smaller 
scales (lA rrigoni Battaia et al.l 120151 : IVerhamme et all 120061 : 
ICantalupo et all 2005lh 

The double peaked Lya profile of companion A, the 
quasar and the bridge (Fig. [5j may be indicative of a sur¬ 
rounding absorbing large neutral gas cloud or h alo, as is 
commonly found around oth er quasars (e.g. iHennawi et al.l 
120061 : IProchaska et al.ll2013al . at similar redshifts). The ab¬ 
sorption in all three components is at the same wavelength 
and goes nearly completely black implying that the absorb¬ 
ing material is separate to the emitting material rather than 
mixed in with it. This absorbing neutral halo has a radius 
> 50 kpc centered on the quasar and hence a mass of at least 
~2 x 10 3 Mq (using the column density estimate above). 
The column density of neutral material in this halo could 



Distance from Quasar /physical kpc 

Figure 6. A comparison of the Lya luminosity along a line 
through the quasar and companion A (taking a cross section along 
the axis of the bridge). The line that seems to fit the data is a 
1/r 2 convolved with a Gaussian of width equal to the PSF (dashed 
line), although there is no physical motivation for this. The fact 
that the surface brightness rises towards A and the quasar sug¬ 
gests that both the active galaxies are sources of the photoioniza¬ 
tion of the bridge. 


be significantly higher than estimated above if in addition to 
the low column density material there are regions of high col¬ 
umn density gas with lower covering fractions. These higher 
density clumps would likely emit Lya photons given that 
they are probably optically thick clouds within the ionizing 
field of the quasar and the companion AGN. However, we do 
not detect a symmetric Lya halo centered on the quasar, as 
opposed to the asymmetric bridge emission, suggesting the 
mass in higher density clumps is significantly less than that 
in the bridge. In terms of Lya absorption, the quasar’s envi¬ 


ronment is consistent with the ty pical environmen t of hig h 
redshift quasars characterised by IProchaska et alJ d2013a ). 
For the bridge and A to show similar absorption they would 
most likely have to be observed through a significant frac¬ 
tion of the absorbing halo, of a comparable column density 
to that affecting the quasar line of sight. 

Although radiative transfer effects can also produce 
double peaked Lya emission in optically th ick gas clouds 
dCantalupo et al' 120051 : [kollmeier et al.luOlOl) . it seems im¬ 
probable that this is the case here as the absorption wave¬ 
length is the same across the entirety of the bridge, quasar 
and companion A whilst the underlying Lya emission profile 
changes significantly. Radiative transfer effects are unlikely 
to produce a double-peaked profile along the line-of-sight 
to the quasar, a point-source, and so the profile here most 
likely arises from multi-phase gas that includes neutral gas 
clouds in front of the quasar. The fact that the absorption 
profile directly in front of the quasar is so similar to that 
across the rest of the system - there are no discontinuities 
in the absorption wavelength - is highly suggestive that the 
quasar line-of-sight is probing the same relative distribution 
of material as the line-of-sight to the bridge. Indeed we see no 
other absorption lines in the quasar spectrum indicating that 
this covering neutral gas must have a low column density, 
consistent with the low column density estimate derived ear¬ 
lier assuming the quasar’s double-peaked profile arises from 
absorbing material rather than radiative transfer effects. 
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The PKS1614+051 system is similar to that of 
SDSSJ0841+3921 and its three companion AGN that are 
also embedded within a n exten ded Lya halo, although at 
slightly lower redshift dHennawi et al . 2mH). However, there 
is no evidence of merging activity or interaction in the 
SDSSJ0841 system and the Lya spectral profiles from the 
quasar and extended nebula are not double-peaked, suggest¬ 
ing there is just one component of (multi-phase) gas present, 
unlike PKS1614+051 that seems to have additional ISM-like 
material inbetween the galaxies. The neutral halo surround¬ 
ing PKS1614+051 may well fluoresce in Lya in a similar 
manner to SDSSJ0841+3921 but at a lower surface bright¬ 
ness than our detection limit. 

Although the picture we put forward here fits with all 
the data it is not necessarily the only interpretation. These 
results are achieved with only 2 hrs of exposure. A signif¬ 
icantly longer observation with MUSE would increase the 
signal-to-noise on the extended Lya emission line allowing 
the spatial and velocity structure of the bridge to be de¬ 
termined in more detail, and potentially detect other emis¬ 
sion lines from the bridge, allowing a plausible model of the 
photoionization source(s) to be put forward. Additionally, 
observations of Oil, OIII and Ha with JWST/NIRSPEC in 
the future would allow the metallicity and ionization param¬ 
eter of the gas to be determined, potentially constraining the 
origin of the gas. 

As noted above, recent work on the gaseous environ¬ 
ments of quasars, in particular the extent of surround¬ 
ing neu tral haloes, has use d background probes (e.g. a 
quasar 


HennawietalJ 120061 : IProchaska &; Hennawil 120091 : 


IProchaska et al]|2013bn to explore the extent of neutral gas 
in the target quasar’s environment. The potential advan¬ 
tages of the approach used here in addition to the use of 
a background probe is that MUSE can allow a search for 
absorption against the extended emission on 10s kpc scales, 
allowing the whole region to be explored rather than just a 
single sight line. The probability of identifying a background 
probe at small projected distance for a quasar decreases as 
(radius)~ 2 so the approach used here can also better sam¬ 
ple the haloes at small radii. Clearly there is potential for 
observing a sample of quasars with extended Lya emission 
with MUSE in order to explore the small s cale quasar envi¬ 
ronmen t. This is particularly important as IProchaska et ahl 
(l2013al l find an anisotropic distribution of N hi absorbers 
around quasars, with more absorbing clouds than expected 
found around the quasar than in its line of sight (the prox¬ 
imity effect), that they attribute to the quasar emitting 
anisotropically - clouds along its line of sight are photo- 
evaporated but those transverse to it are exposed to less ra¬ 
diation and survive. We see the Lya absorption in PKS1614 
across the whole of the quasar and companion system sug¬ 
gesting that this quasar is radiating relatively isotropically. 
In addition, given that we see the bridge and both AGNs, 
and the bridge must lie mainly in the plane of the sky or 
its true extent becomes improbably large, the opening angle 
for AGN emission in both A and the quasar must be large 
and any obscuring torus small. 

The profile of the Lya absorption is symmetric around 
the quasars systemic velocity suggesting it is not associated 
with any inflow or outflow of material. If material was in¬ 
flowing or out-flowing it would need to be continually replen¬ 
ished or destroyed, or be a short lived phenomenon, as the 


amount of absorber material is relatively small (> 2 x 10 3 
M@). Instead we postulate this neutral material sits spher¬ 
ically around the quasar within its dark matter halo and 
(some of) its large velocity width comes from the fact that 
these clouds are moving within the gravitational potential 
of the halo. It is clear that MUSE, with its large IFU cov¬ 
erage, provides an excellent dataset to probe the quasar’s 
environment over an extended region, including providing 
important constraints on the properties of its halo gas. 


5 CONCLUSIONS 

We find that the z = 3.215 quasar PKS1614+051, with its 
companion active galaxy and connecting bridge of Lya emit¬ 
ting gas, is surrounded by three other Lya emitting galaxies, 
two of which are very close to or within the bridge of mate¬ 
rial. This overdensity of four galaxies, all within 300 km s _1 
of each other and within the virial radius of the quasar, will 
likely merge to form a massive elliptical galaxy by 2 ~ 1, 
perhaps residing within a group or cluster of galaxies. 

The 50 kpc long bridge of Lya emitting gas stretches be¬ 
tween the quasar and companion galaxy A and also extends 
towards the two other nearby companion galaxies, suggest¬ 
ing that the gas was pulled out from the galaxies during a 
gravitational interaction between them and the quasar host. 
The emission line bridge is likely to be predominantly pho- 
toionized by the quasar and its AGN companion as its sur¬ 
face brightness rises towards them both, although the lack of 
CIV or Hell in the bridge suggests a relatively low ionization 
parameter. The quasar’s Lya line is strongly double-peaked, 
likely due to absorption by a neutral halo of gas at the sys¬ 
temic velocity of the quasar extending at least 50 kpc in 
radius with a mass of at least ~2 x 10 3 Mq. The Lya bridge 
is probably a separate phase to the much larger neutral halo 
surrounding the quasar, which it is embedded in, as the Lya 
absorption profile goes nearly completely black indicating 
virtually all of the Lya photons at that wavelength have 
been scattered out of the line of sight. In addition, we must 
be looking at the system through a significant path length of 
neutral material to see similar absorption across the entire 
system. This picture of the quasar, bridge and companion 
galaxies embedded in a large neutral halo is consistent with 
the absorption systems commonly found ar ound quasars us¬ 
ing a b ackground quasar as a probe by e.g. IProchaska et al.l 
(|2013aT ). The fact that we see the bridge and both of the 
AGNs, and that we see neutral material in both the line of 
sight to the quasar and transverse to it, suggests the quasar 
is emitting isotropically and any obscuring torus obstructs 
only a small solid angle as seen from the nucleus. 

The large field of the MUSE IFU allows the exploration 
of the quasar environments on a wide range of scales. The 
clustering of galaxies at the quasar redshift can be stud¬ 
ied simultaneously with the small-scale exploration of the 
quasar’s gaseous environment using Lya absorption of faint 
sources in the immediate environment. This use of self¬ 
absorption, rather than using a rare background quasar to 
probe a single line of sight, is a potentially powerful tech¬ 
nique for exploring the close environments of quasars, given 
the prevalence of extended Lya emission. 

This work demonstrates the effectiveness and efficiency 
of MUSE for elucidating the observational details of galaxy 
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interactions and growth at high redshift. While a single 
pointing can explore the early growth of clustering, mini¬ 
mal tiling of several MUSE pointings should ideally sample 
the larger-scale protocluster environment at these redshifts. 
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